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ABSTRACT 

We expand our Bayesian Monte Carlo method for analyzing the light curves of 
gravitationally lensed quasars to simultaneously estimate time delays and quasar struc- 
ture including their mutual uncertainties. We apply the method to HE1104-1805 and 
QJ0158-4325, two doubly-imaged quasars with microlensing and intrinsic variability on 
comparable time scales. For HE1 104-1805 the resulting time delay of AtAB = iA — ^B = 
162.2jjg days and accretion disk size estimate of log(r s /cm) = 15.7i ' 5 at 0.2/iin in the 
rest frame are consistent with earlier estimates but suggest that existing methods for 
estimating time delays in the presence of microlensing underestimate the uncertainties. 
We are unable to measure a time delay for QJ0158-4325, but the accretion disk size is 
log(r s /cm) = 14.9 =L 0.3 at 0.3/xm in the rest frame. 

Subject headings: cosmology:observations — accretion, accretion disks — dark matter 
— gravitational lensing — quasars: general 
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Introduction 



Variability in lensed quasar images comes from two very different sources. Changes in the 
quasar's intrinsic luminosity are observable as correlated variability between images, while mi- 
crolensing by the stars in the lens galaxy produces uncorrelated variability. Measurements of the 
time del ays between th e lensed ima ges from the correlated varia bility can be used to study cosmol- 



ogy (e.g. lRefsdallll964l and recently 



in the lens galaxy (e.g. iKochanek et al 



Saha et al 
2006; 



2006 



Poindexter et al 



Ogu rill2007r) or the distribution of da rk matter 



2007 



Vuissoz et alJl2007l >. The mi- 



crolensing variability can be used to study the structure of the quas ar, the masses of the stars in the 
lens galaxy, and th e stellar mass fraction near the lensed images (jSchechter k, Wambsganss 2002; 
Wambsganss! 120061 ) . It is now possible to use micro lensing to measure the correlation of accretion 
disk size with b lack hole mass ([Morgan et al.l 120071). th e wavelength dependence of the size of the 
accretion disk ( Poindexter. Morgan KochanekT 2007 ) or the differing sizes of the thermal and 
non-thermal X-ray emission regions (jPoolev et alj|2007l : iDai et ajJ|2007J ). 



The challenge is that most lensed quasars exhibit both intrinsic and microlensing variability. 
To measure a time delay, one must successfully model and remove the microlensing variability 
such that only intrinsic variability remains. If the microlensing variability has a suffi ciently low 
amplitude or long timescale, it can be igno red (e.g. PG1 1 15-080. ISchechter et alJll997l ). but this is 
a dangerous assumption for many systems. Eigenbrod et al. ( 20051 ) found that for an 80 day delay, 
adding microlensing perturbations with an amplitude of 5% (10%) to a light curve increased the 
uncertainty in th e time delay by a factor of 2 (6). Existing time delay analyses for lenses with 
microlensing (e.g. iParaficz et al.ll2006l ; IKochanek et al.ll2006l ; iPoindexter et al.ll2007l ) depend on the 
intrinsic and microlensing variability having different time scales. These analyses also require that 
the microlensing variability can be modeled by a simple polynomial function. This approach will 
clearly fail if the two sources of variability have similar time scales or if the microlensing variability 
cannot be easily parameterized. 

In this paper, we present a new technique for simultaneously estimating the time delay and 
structure of lensed quasars that exhibit strong microlensing. In essence, we assume a range of 
time delays and then determine the likelihood of the impli ed microlensing variab ility using the 
Bayesian Monte Carlo method of Kochanek (2004, see also IKochanek et al.ll2007l ). This allows 
us to estimate the time delays and the quasar structural parameters simultaneously and include 
the effects of both phenomena on the parameter uncertainties. We apply the method to the two 
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doubly-imaged lenses HE1104-1805 (IWisotzki et alJll993T) and QJ0158-4325 (|Morgan et alJll99a ). 
While HE1104-1805 has a well-measured time delay (jPoindexter et al.ll2007l ). the amplitude of the 
microlensing (~ 0.05 mag yr" 1 over the past decade) and the fact that it exhibits variability on the 
6 month scale of the time delay suggest that it is close to the limit where microlensing polynomial 
fitting methods (jBurud et al.ll200ll ; iKochanek et al.l 12000 ) will break down. QJ0158-4325 clearly 
shows both correlated and uncorrelated variability, but the polynomial methods cannot reliably 
produce a time delay estimate. We describe the data and our models in $2j our new approach in $3] 
and the application to the two systems in §H In § El we discuss the results and their limitations. We 
assume a flat = 0.3, Aq = 0.7, Hq = 70 km s _1 Mpc~ x cosmology and that the lens redshift of 
QJ0158-4325 is Z\ = 0.5. Reasonable changes in this assumed redshift have negligible consequences 
for the results. 



2. Observations and Models 



We created photometric models for the two systems from the WFPC2 and NICMOS V- 
(F555W), I- (F814W) and ff-band (F160W) observations of the two systems by the Cf A- Arizona 
Space Telescope Lens Survey (CASTLES) following the methods of lLehar et al.l (120001 ). The quasars 
are modeled as point sources and each lens galaxy as a de Vaucouleurs profile. We chose the de 
Vaucouleurs profile above other models since it provided the best fit. Table Q] summarize s the fits 
we use here, where the HE1 104-1805 model is updated from t hat in iLehar et al.l (|2000l ) using a 
deeper if-band image obtained to study the quasar host galaxy (|Yoo et al.ll2006l ). 



For each system we created a sequence of ten lens models using the lensmodel software pa ckage 
(jKeetonl l200ll ). Each model is the sum of concentric NFW (jNavarro. Frenk fe White! Il996l ) and 
de Vaucouleurs components, where the NFW component simulates the dark matter halo and the 
de Vaucouleurs component represents the galaxy's stellar content. We parameterize the model 
sequence by /m/l> t ne mass of the stellar component relative to its mass in a constant mass- 
to-light ratio model with no contribution from the NFW halo. We generated model sequences 



covering the range 0.1 < Jm/l — 1-0. With their time delay measurement, iPoindexter et al 



Table 1. HST Astrometry and Photometry of QJ0158-4325 and HE1 104-1805 



Lens 


Component 


Astrometry 




Photometry 








ARA 


ADec 


H=F160W 


I=F814W 


V=F555W 


QJ0158-4325 


A 


= 


= 


16.47 ±0.03 


17.81 ±0.04 


18.10 ±0.13 




B 


-l'/156 ± 0'.'003 


-0'.'398 ± 0'.'003 


17.27 ±0.03 


18.62 ±0.11 


18.91 ±0.17 




G 


-0'.'780±0'.'016 


-Of.' 234 ± 0'.'006 


16.67 ±0.13 


18.91 ±0.06 


20.36 ±0.18 


HE1104-1805 


A 


= 


= 


15.91 ±0.01 


16.40 ± 0.03 


16.92 ±0.06 




B 


+2'/901 ± 0'.'003 


-1'.'332±0'.'003 


17.35 ±0.03 


17.95 ±0.04 


18.70 ±0.08 




G 


+0'/965 ± 0'.'003 


-0'.'500 ± 0'.'003 


17.52 ±0.09 


20.01 ±0.10 


23.26 ±0.27 
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Fig. 1. — Top Panel: QJ0158-4325 R-band light curves for images A (squares) and B (triangles). 
The curves are a polynomial fit to guide the eye. Points with error bars only (no symbols) were not 
used in the analysis. Lower Panel: Example of a good fit to the implied A-B microlensing signal 
for a trial time delay of Ai^B = — 20 days. 
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( 20071 ) constrained the stellar mass fraction of HE1104-1805 to f M/L = 0.30±g;g|, aut we chose 
not to apply these limits to /m/l f° r our present calculations. From these models for the mass 
distribution we extract the convergence (k), shear (7) and stellar surface density (k*) for each image 
and then generate realizations of the micr olensing magnification patterns at the locati on of each 
i mage using a variant of the ray-shooting ([Schneider et al.l [1992]) method described in iKochanek 
(200J). We assume a stellar mass function of dN(M)/dM oc M~ 1 ,3 wit h a dynamic range of 50 



which approximates the Galactic stellar mass function of iGouldl (12000 ) . We present the results 
either making no assumption about the mean mass (M) of the stars or by applying a prior that 
it lies in the range 0.1 M < (M) < 1.0 Mq. We used 4096 2 magnification patterns with an outer 
scale of 20Re- We used a prior for the relative motions of th e observer, lens ga laxy, lens galaxy stars 
and the source based on the projection of the CMB dipole (jKogut et al.lll993l ) for the observer, the 
stellar velocity dispersion of the lens set by its Einstein rad ius, and rms peculiar velocities for the 



lens and source of a p = 235/(1 + z) k m s 1 ( Kochanek 2004). We modeled the continuum emission 



source as a face-on thin accretion disk (jShakura &: Sunvaevlll973l ) with the surface brightness profile 



I(R) oc {exp [{R/r s f 4 } - l} * . 



(1) 



The scale radius r s is the point where the disk temperature matches the rest frame wavelength of 
our monitoring band, kT = hc(l + z s )/X f JS - For comparisons t o othe r disk models, the half-light 
radius of R1/2 = 2.44r s should be used, since iMortonson et al.l (|2005l ) have shown that half-light 
radii estimated from microlensing depend little on the assumed surface brightness profile of the 
disk. 

We monitored QJ01 58-4325 in the fl-ba nd using the SMARTS 1.3m telescope with the ANDICAM 
optical/infrared camera (|DePov et alJhooal ff and using the 1.2m Euler Swiss Telescope CIS £1 part 



of the the CQS MOGRAIljl 



can be found in Kochanek et al. 



jaro 



, e ct. A full description of our monitoring data reduction technique 
(2006), but we provide a brief summary here. We model the PSF 
of each quasar image using three nested, elliptical Gaussian components, keeping the relative as- 
trometry fixed for all epochs. We use relative photometry, comparing the flux of each image to the 
flux of several reference stars in each frame. The lens galaxy flux is determined by optimizing its 
flux in observations with good seeing, and the light curves are then measured with the galaxy flux 
fixed to this optimal value. We eliminated all data points taken at seeing conditions worse than 
l'/7. We dropped 3 epochs that satisfied the seeing conditions but on which the sky was very bright. 
On these nights, the flux measurements from the dimmer image B are clearly contaminated by sky 
flux, so we report these measurements but chose not to include them in our microlensing analysis. 
The monitoring data are presented in Table [21 and the light curves are d isplayed in Figure H Fo r 
HE1104-1805, we use the composite R- and V-band light curve data from lPoindexter et al.l (|2007l ). 



2 http: //www. astronomy, ohio-state.edu/ ANDICAM/ 
3 http://www. cosmograil.org/ 
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Fig. 2. — Probability distribution for the time delay of HE1104-1805, where AtAB =tA~ ^B- The 
solid curve has no prior on the microlens mass scale while the dotted curve assumes a uniform prior 
on the mass over the range 0.1 Mq < (M) < l.OMp. Previou s measurements of the time delay 
are also plotted. The dependence of the iPoindexter et al.l (120071 ) results on the statistical test used 
demonstrates the limitations of polynomial fitting methods. 
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Fig. 3. — Probability distributions for the source size r s of HE1104-1805. The dashed curve shows 
the estimate for r s with a prior of 0.1 < (M) < 1.0 M Q on the mean mass of the microlenses. 
The vertical line shows the Schwarzschild radius Rbh = 2GMbh I ' ( ? of t he black hole, where the 
black hole mass M B h = 2.37 x 1O 9 M was estimated by iPeng et al.1 ()2006l ) using the C IV emission 
line width. The last stable orbit for a Schwarzschild black hole is at 2>Rbh- 
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3. Joint Monte Carlo Analysis 

For our analysis we must generate light curve pairs with arbitrary delays. We always carry out 
the shifts on the less variable light curve (image A for both HE1104-1805 and QJ0158-4325). There 
are two issues for generating the light curves. First, we must define the algorithms for interpolating 
and extrapolating a light curve, and second we must decide how many extrapolated points should 
be used. When data points are shifted in the middle of an observing season, we estimate the flux at 
the shifted time using linear interpolation between the nearest bracketing data points. When shifted 
points lie in an inter-season gap or beyond either end of the observed lightcurve, we estimate the 
flux using extrapolation based on a linear fit to the five nearest data points. Data points requiring 
extrapolation for periods longer than seven days are discarded. We assign new uncertainties that 
combine the photometric errors with the uncertainties due to the temporal distance of the new 
point from existing data point s , whe re we model the source variability using the quasar structure 



function of IVanden Berk et al.l (j2004l ) 



For a standard analysis, we wanted to use light curves with the same number of data points 
for each time delay and no points extrapolated by more than seven days. We limit the use of 
extrapolated points because they lead to a delay-dependent change in the statistical weights caused 
by the steadily growing uncertainties from the structure function. Given the delay range we wanted 
to test, we first found the limiting delay defined by the delay which yielded the minimum number 
of usable data points given our seven day extrapolation limit. We then restricted all of the trial 
light curves to use only the epochs permissible at the limiting delay. For the very long delays of 
HE1104-1805, the limiting delay is not the longest delay, because the longest delays shift curves 
completely through the inter-season gaps. This forced us to restrict the HE1104-1805 analysis to 
the epochs permitted by both the limiting delay and the longest delay, since some permissible times 
at the end of the limiting delay lightcurve are beyond the extrapolation limit in light curves with 
longer delays. We also experimented with simply allowing unlimited extrapolations. 



We then analyzed the light curves using the Bayesian Monte Carlo method of lKochanekl (|2004l ) . 
In essence, we randomly select a time delay At, a lens model and a disk model (size), generate a 
microlensing light curve and fit it to the microlensing light curve implied by the observed data and 
the selected time delay. This gives us a x 2 -statistic for the goodness of fit for the trial ^{p^At) 
given the model parameters for the microlensing p = (/m/l, r s-> velocities, masses etc.) and the 
time delay. Figure Q] shows an example of a good trial lightcurve fit to QJ 0158-4325 for a delay of 
AtAB = —20 days. In essence, the probability of time delay At is the Bayesian integral 

P{At\D) oc J P{D\p,At)P(p)P(At)dp (2) 

where P(D\p, At) is the proba bility of fitting the data in a particular trial, P{p) sets the priors on 



the microlensing variables (see iKochanek! 12004! : iKochanek et al.ll2007l ) and P{At) is the (uniform) 



prior on the time delay. The total probability is then normalized so that J P(At\D)dAt = 1. 
We evaluated the integral as a Monte Carlo sum over the trial light curves, where we created 
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Fig. 4. — Probability distribution for the time delay AtAB = tA — *b i n QJ0158-4325, where 
we restrict the extrapolation of the light curves to be less than 7 days. The solid curve assumes 
no prior on the microlens mass and the dashed curve assumes a uniform prior over the range 
O.1M < (M) < l.OM . 
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Fig. 5. — Same as Fig. [4J but the trial light curves were generated with no limits on extrapolation. 
Note the change in the relative probability of positive and negative delays. 
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4 independent sets of magnification patterns for each of the ten macroscopic mass models and 
generated 4 x 10 6 trial light curves for each magnification pattern set. 



4. Results 



Poindexter et al.1 (120071 ) recently estimated a time delay for HE1104-1805 of AtAB = tA — tB = 
152.2 + 3'q day s, in t he se nse that image A lags ima ge B, improving on the earlie r estimates by 
Ofek fc Maozl ||2003h and IWvrzvkowski et al.1 (|20o3 ). The IPoindexter et al.l (j2007h analysis used 



the 



of Schechter et al 



Kochanek et al,l ( 20061 ) polynomial method on li ght curves which combined the published data 



|2003h . bfek fc Mao3 feoO$ ) and IWvrzvkowski et all d2003h with new i?-band 
monitoring data. In this polynomial method, the source and microlensing variability are modeled 
as a set of Legendre polynomials that are then fit to the light curves. Ambiguities arise because the 



value of the delay depends weakly on the parameterization of the microlensing. IPoindexter et al 



(|2007l ) used a Bayesian weighting scheme for the different polynomial orders, but obtained 157.2 ± 
2.6 days if they used the F-test to select among the different orders rather than a Bayesian weighting. 
The advantage of our present approach is that it uses a physical model for the microlensing rather 
than a polynomial parameterization of it. 

We applied our joint Monte Carlo analysis technique to HE1104-1805 over a time interval 
range of 125 days < AtAB < 200 days with a sampling in 1.5 day intervals and no extrapolation of 
the light curves past 7 days. Figs. [2] and [3] show the resulting probability distributions for the time 
delay and the disk size. We find a ti me delay of AtAB =_tA — ts = 162.2^'g days (la) that is in 



marginal agreement with the formal IPoindexter et al.1 (120071 ) result but in better agreement with 
the F-test selection of the best polynomial model than with the Bayesian result. Assuming Hq = 
70 kms -1 Mpc -1 and an inclination angle of cos(i) = 1/2, the disk size estimate of log(r, 9 /cm ) = 



15.7 + Q'g at 0.2 um in the rest frame is lit tle changed from the estimates in Morgan et al.l (12007T) and 



Poindexter. Morgan Kochanekl (|2007l ) which held the delay fixed to the IPoindexter et al.l (120071 ) 
value. The uncertainties in the size are a factor of 1.6 larger in the present analysis. 

In models of QJ0158-4325, we expect image A to lead image B by 2-3 weeks depending on the 
mass distribution and the actual lens redshift (we assumed zi = 0.5). For completeness, we tested 
a full range of negative and positive delays —50 days < AtAB < 50 days at 2 day intervals. We 
generated the positive and negative delay light curves separately in order to minimize the number 
of points lost due to the 7 day extrapolation limit at the price of making the data used for positive 
and negative delays somewhat different. For comparison, we also tried using all points with no 
extrapolation limits. 

Figs. H] -[6] show the results for the time delays and the source size. We have clearly failed to 
measure a time delay, and positive delays AtAB = tA~tB seem to be favored over negative delays, in 
direct disagreement with the predictions of the lens model. The relative likelihoods of positive and 
negative delays depend on the detailed treatment of the light curves, with the probability of positive 
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Fig. 6. — Probability distributions for source size r s of QJ0158-4325. The dashed (dotted) curve 
shows the estimate for r s using the set of negative (positive) trial time delays. The solid curve 
is the r s estimate using all trial time delays. The vertical line shows the Schwarzschild radius 
Rrh = 2GMrh /c 2 o f the black hole, where the black hole mass Mbh = 1-6 x 1O 8 M0 was estimated 



bv lPeng et al.l (|2006l ) using the Mg II emission line width. The last stable orbit for a Schwarzschild 
black hole is at 3Rbh- 
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delays being lower when we use all extrapolated points rather than restricting them to 7 days or 
less. In this case, the microlensing simply overwhelms the intrinsic variability. We expect additional 
monitoring data to continue to tighten the time delay probability distribution, but a successful delay 
measurement may not be possible for many seasons. We succeed, however, in estimating a size for 
the quasar despite the uncertainties in the time delay, finding a size of log(r s /cm) = 14.9 ± 0.3 
at 0.3 fim in the rest frame (again, assuming Hq = 70 kms -1 Mpc -1 and cos(z) = 1/2). The size 
estimate changes little compared to the uncertainties if we limit the analysis to either positive or 
negative delays. The microlensing amplitudes are large enough that interpreting varying amounts 
of the lower amplitude intrinsic variability as microlensing does not change the statistics of the 
microlensing enough to significantly affect the size estimate. One additional uncertainty in this 
result is that lens redshift of QJ0158-4325 is unknown. We experimented with running the Monte 
Carlo simulation at a range of lens redshifts (0.1 < z\ < 0.9), and we found that the resulting shifts 
in the r s estimates were negligible relative to the size of the existing uncertainties. 



5. Discussion and Conclusions 



Peng et al.1 (1200fih used the width of the C IV (A1549A) emission line to estimate the black 
hole mass M BH = 2.37 x 1O 9 M in HE1104-1805 and the width of Mg II (A2798A) emission line 
to estimate the black hole mass Mbh = 1-6 x 10 8 M Q in QJO 158- 4325. Using these b lack hole 
masses, the quasar accretion disk size - black hole mass relation of Morgan et al.1 ( 2007 ) predicts 
source sizes at 2500A of log(r s /cm) = 15.9 ± 0.2 for HE1104-1805 and log(r s /cm) = 15.2 ± 0.2 for 
QJ0158-4325. If we scale our current disk size measurements to 2500A using the R\ oc A 4 / 3 scaling 
of thin disk theory and assume an inclination angle cosi = 1/2, we find log(r s /cm) = 15.91q' 4 for 
HE11 04-1805 and log(r, s /c m) = 14.8 ± 0.3 for QJ0158-4325, fully consistent with the predictions 
of the Morgan et al.l (120071 ) accretion disk size - black hole mass relation. 



The mixing of intrinsic and microlensing yariability in lensed qu asar light curves can be a 
serious problem for estimating time delays (e.g. lEigenbrod et al.ll2005l ) and previous microlensing 
analyses have been restricted to lenses with known time delays. In HE1104-1805, which must 
be close to the limits of measuring time delays in the presence of microlensing, we confirm that 
the approach of fitting polynomial models for the microlensing works reasonably well. However, 
the dependence of the delay on the assumed model was a warning sign that the formal e rrors 
on the delays were likely to be underestimates, as was recognized by iPoindexter et al.1 (|2007l ). In 
our new, non-parametric microlensing analysis of HE1104-1805 we find a modestly longer delay of 
162.2^'g days that quantifies those concerns. Estimates of the quasar accretion disk size are little 
affected by these small shifts in the time delay. In QJ0158-4325, the microlensing amplitude is 
larger relative to the intrinsic variability, and traditional methods for determining delays fail. Our 
new method also fails to measure a delay, but it does allow us to measure the size of the quasar 
accretion disk despite the uncertainties in the time delay. 
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Table 2. QJ0158-4325 Light curves 



HJD X 2 /N do f QSO A (mags) QSO B (mags) (Stars) Source 
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blvlAKl b 


3677. 


bbl 


o 



A 1 

.41 


1 


noo 
.922 


± 0. 


nn^7 

.007 


2 


orr 

.000 


± 0. 


,009 


0. 


1 O 1 

131 


± 0. 


,002 


T7TTT IT O 


3678. 


.b44 


-i 

1. 


.28 


1 


non 

,9o9 


± 


,009 


2 


onn 

.291) 


± 0. 


Oil 


— 0. 


.004 


± 0. 


nno 

.003 


O A /T A T ) C 

blvlAKl b 


0081 


f A O 

.643 


u. 


o o 
,55 


1 


nro 

.958 


i n 

± 


nnn 

.009 


2 


.345 


i n 

± 0. 


m o 
,1)12 


— 0. 


noo 

033 


i n 

± 


nno 

.003 


O A H A TDT^O 

blvlAKl b 


0050 


f f n 

.669 


r 

5. 


1 A 

,14 


1 


.914 


i n 

± 


nnf 

.006 


2 


A A<~1 

.442 


i n 

± 0. 


nno 

,008 


n 

0. 


nTn 

.079 


i n 

± o. 


nno 

.003 




Of O O 

0050. 


f O O 

.633 


U. 


,50 


1 


.944 


i n 

± 


n 1 1 

.011 


2 


O TO 

.373 


i n 

± 0. 


,015 


— 0. 


nf a 

.064 


i n 

± 


nno 

.003 


OA /T A T^T^O 

SMAK1 S 


Of o o 

0055. 


f Tf 

.676 


4. 


.26 


1 


.910 


i n 

± 


nnf 

,006 


2 


ATI 

.451 


i n 

± 0. 


nno 

,008 


0. 


nf n 

.069 


i n 

± 


nno 

,003 




of no 

3692 


f O 1 

.631 


2. 


— o 

,72 


1 


noo 

.923 


i n 

± 


nnf 

.006 


2 


/I Of 

.436 


i n 

± 0. 


nno 

,008 


0. 


nn /i 

.094 


i n 

± 


nno 

.003 




of no 

3693 


tr f n 

.560 


2. 


1 n 

.10 


1 


n /i i 

.941 


i n 

± 


nno 

.008 


2 


/i nn 

.409 


i n 

± 0. 


m n 

,010 


— 0. 


.012 


i n 

± 0. 


nno 

.003 


HIT* DTP 

SMAK1 S 


3694 


.631 


2. 


.43 


1 


.946 


± 


.007 


2 


.455 


±0. 


,009 


0. 


113 


± 


.002 


EULER 


3696 


.623 


5. 


.99 


1 


.923 


±0 


,006 


2 


.503 


±0. 


.008 


0. 


.135 


±0. 


.002 


EULER 


3700 


.648 


2. 


.38 


1 


.931 


±0 


.006 


2 


.452 


±0. 


.009 


0. 


.122 


±0. 


.002 


EULER 


3701 


.605 


1 


.01 


1 


.951 


±0. 


.010 


2 


.400 


±0. 


.013 


-0. 


021 


±0. 


,003 


SMARTS 


3705 


.652 


1, 


.12 


1 


.952 


±0 


.009 


2 


.389 


±0. 


.011 


-0. 


.004 


±0. 


.003 


SMARTS 


3707 


.706 


2, 


.82 


1 


.951 


±0 


.007 


2 


.450 


±0. 


.010 


0. 


,122 


±0. 


,002 


EULER 


3710. 


.605 


0. 


.99 


1 


.957 


±0 


.010 


2 


.424 


±0. 


.014 


-0. 


012 


±0. 


.003 


SMARTS 


3715. 


693 


2, 


.21 


1 


.967 


±0. 


.006 


2 


.452 


±0. 


.007 


0. 


,105 


±0. 


,002 


EULER 


3717. 


.573 


0. 


.68 


1 


.953 


±0 


.010 


2 


.459 


±0. 


.015 


-0. 


061 


±0. 


.003 


SMARTS 


3720. 


.644 


3. 


.80 


1 


.949 


±0. 


.007 


2 


.480 


±0. 


.009 


0. 


.062 


±0. 


,003 


EULER 


3732 


.631 


4. 


.36 


1 


.999 


±0 


.007 


2 


.516 


±0. 


.010 


0. 


.121 


±0. 


.002 


EULER 


3735 


.616 


4, 


.39 


1 


.985 


±0. 


.006 


2 


.541 


±0. 


.009 


0. 


.126 


±0. 


,002 


EULER 


3747. 


,614 


3. 


.03 


1 


.990 


±0 


.006 


2 


.576 


±0. 


.008 


0. 


079 


±0. 


.003 


EULER 
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Table 2 — Continued 



HJD X 2 /N do f QSO A (mags) QSO B (mags) (Stars) Source 



ott t 

3757. 


.606 


1 

1 


no 

.92 


1 

1 


nno 

.998 


i n 

± 


nno 

,008 


2 


r n 1 

,591 


i n 

± 0. 


m o 

.013 


0. 


no o 

,088 


i n 

± 0. 


nno 

.003 


T7>TTT TA T") 


3764. 


.534 


n 




on 
,81) 


1 


.994 


± 0. 


m n 

.010 


2. 


.579 


i n 
± 0. 


ni £ 

.015 


n 
— 


.033 


i n 

± 0. 


nno 

.003 


MYlARl S 


3765. 


.570 


r- 



no 

.11 


2 


nn^7 

.007 


± 0. 


.006 


2. 


.575 


i n 
± 0. 


.007 


n 



1 on 

.130 


i n 

± 0. 


nno 

.002 


T^TTT T7 D 


O '7*71 

3771. 


.604 


2 


O 1 

.31 


2 


nno 

.003 


i n 

± 0. 


.006 


2. 


r n o 

,598 


i n 

± 0. 


nnn 

.009 


n 




nn 1 

.091 


i n 

± 0. 


nno 

.003 


T^TTT T7 D 


3782. 


.557 


1 


c n 

,59 


1 


no a 

.984 


i n 

± 0. 


nno 

.008 


2. 


.565 


i n 

± 0. 


m o 

.013 





non 

.032 


i n 

± 0. 


nno 

.003 


T^TTT T7 T } 
IIjULIIjK 


3787. 


.531 


2 


.82 


2 


.027 


± 0. 


.006 


2. 


.549 


± 0. 


.009 





.098 


± 0. 


.003 


TATTT T - " T ~) 

hULhR 


3800. 


n a 

.519 


a 

2 


no 

.23 


a 

2 


n i 1 

.011 


i n 

± 


n i n 

.010 


A 

2 


.576 


i n 

± 0. 


A-i r 

.015 


A 

0. 


non 

,032 


i n 

± 0. 


nno 

.003 


TATTT T - " T~) 

hULhR 


3889. 


non 

.922 


2 


.68 


1 


non 

.980 


i n 

± 0. 


nno 

.008 


A 

2 


r r 1 

.551 


i n 

± 0. 


m o 

.013 


n 

0. 


AAf7 

.097 


i n 

± 0. 


nno 

.003 


T7>TTT T7 , T~> 


3908. 


A1 A 

910 


2 


A A 

.44 


1 


n r o 

.953 


i n 

± 


AAA 

.006 


A 

2. 


r /i o 

,543 


i n 

± 0. 


nno 

,008 


n 

u. 


1 no 

.128 


i n 

± 0. 


nnn 

.002 


T7>TTT TT 1 T~) 

hULhR 


3913 


.839 


2 


n ^ 

.94 


1 


no*7 

.987 


i n 

± 0. 


nn*7 

.007 


2. 


rrr 

.555 


i n 

± 0. 


nnn 

.009 


0. 


1 no 

.102 


i n 

± 0. 


nno 

.003 


T7>TTT T7 , T~> 

hULhR 


3919. 


oon 

.880 


5 


n i 

.01 


-I 

1 


n r 

.956 


i n 

± 


nnn 

,009 


A 

2 


ror 

,585 


i n 

± 0. 


m a 

.014 


n 

0. 


nnn 

,029 


i n 

± 0. 


nno 

.003 


T7>TTT T7 T~) 

hULhR 


3930. 


or/ 1 

85b 





a n 

.49 


1 


no o 

.983 


i n 

± 0. 


nn 1 

.021 


A 

2 


r 77 1 

.577 


i n 

± 0. 


AOA 

.036 


A 

— 0. 


A1 A 

.016 


i n 

± 0. 


nn a 

.004 


S1V1AK1 S 


3932 


.920 


1 


1 n 

.40 


1 


n7o 

.973 


i n 

± 


/VI 1 

.011 


A 

2. 


,611 


i n 

± 0. 


m n 

,019 


n 

0. 


.007 


i n 

± 0. 


nno 

,003 


TT^TTT T7 T") 


3950 


oon 

.832 


2 


.46 


1 


nno 

.928 


i n 

± 0. 


AAA 

.006 


A 

2. 


r on 

.539 


i n 

± 0. 


nno 

.008 


n 

0. 


i /i n 

.140 


i n 

± 0. 


nnn 

.002 




3960 


O 1 7 

.817 


1 


no 

.23 


1 


nn*7 

.997 


i n 

± 0. 


n i n 

.010 


A 

2. 


.552 


i n 

± 0. 


m a 

.014 


A 

— 


.051 


i n 

± 0. 


nno 

.003 


S1V1AK1 S 


3961 


nnn 

.922 


6 


,26 


1 


n A a 

.944 


i n 

± 0. 


nn/> 

.006 


A 

2. 


.594 


i n 

± 0. 


nnn 

.009 


n 

0. 


nno 

.093 


i n 

± 0. 


nno 

.003 




3967 


.838 





,85 


1 


.957 


± 0. 


.010 


2 


,546 


± 0. 


.014 


— 


,026 


± 0. 


.003 


OT\ if A T~lT~iO 

SMAK1 S 


3974. 


.783 


1 


,11 


1 


.984 


± 0. 


.010 


2 


.550 


± 0. 


.014 


— 0. 


.030 


± 0. 


.003 


OT\ T A T~*T"iO 

SMARTS 


3995 


.765 





.88 


1 


A A 7 

.927 


i n 

± 0. 


n i 1 

.011 


A 

2. 


.511 


i n 

± 0. 


A1 A 

.016 


— 0. 


no A 

.034 


i n 

± 0. 


nno 

.003 


OT\ T A T"»T"*0 

SJV1AK1S 


< a a a 

4002. 


7nn 

.702 


1 


no 

,03 


1 


no 1 

.931 


i n 

± 0. 


n i o 

.012 


A 

2. 


Af\ A 

,494 


i n 

± 0. 


n i o 

.018 


— 


n /i o 

,043 


i n 

± 0. 


nno 

.003 


OT\ T A TWO 

SMAK1S 


4007 


.696 





,70 


1 


.905 


±0. 


,010 


2 


,479 


±0. 


.015 


-0 


,032 


±0. 


.003 


SMARTS 


4030 


.632 





,82 


1 


.876 


±0. 


.010 


2 


.372 


±0. 


.013 


-0. 


.020 


±0. 


.003 


SMARTS 


4037 


.600 





,92 


1 


.886 


±0. 


.009 


2. 


.390 


±0. 


.012 


-0. 


.025 


±0. 


.003 


SMARTS 


4043. 


.600 





.48 


1 


.861 


±0. 


.015 


2. 


,397 


±0. 


.023 


-0. 


,063 


±0. 


.004 


SMARTS 


4050. 


.620 


1 


.34 


1 


.888 


±0. 


.008 


2. 


.393 


±0. 


.011 


-0. 


.022 


±0. 


.003 


SMARTS 


4062. 


.536 





.47 


1 


.891 


±0. 


.051 


2. 


.432 


±0. 


,084 


0. 


,109 


±0. 


.005 


SMARTS 


4064. 


.646 





.84 


1 


.881 


±0. 


.009 


2. 


.434 


±0. 


.012 


-0. 


.011 


±0. 


.003 


SMARTS 


4069. 


.580 





.82 


1 


.872 


±0. 


.009 


2. 


.419 


±0. 


.013 


-0. 


,049 


±0. 


.003 


SMARTS 


4083. 


.542 





.69 


1 


.895 


±0. 


.010 


2 


.527 


±0. 


.015 


-0. 


.037 


±0. 


.003 


SMARTS 


4090. 


621 


2 


.29 


1 


.865 


±0. 


.010 


2. 


.541 


±0. 


,014 


-0. 


,037 


±0. 


.003 


SMARTS 


4097. 


.556 





.40 


1 


.882 


±0. 


.015 


2. 


.554 


±0. 


.027 


-0. 


.059 


±0. 


.004 


SMARTS 


4111. 


.601 


1 


.23 


1 


.880 


±0. 


.011 


2 


,622 


±0. 


.019 


-0. 


,033 


±0. 


.003 


SMARTS 


4114. 


.561 


1 


,08 


1 


.842 


±0. 


.013 


2. 


.664 


±0. 


.027 


-0. 


.037 


±0. 


.003 


SMARTS 
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Table 2 — Continued 



HJD 


X 2 /N dof 


QSO A (mags) 


QSO B (mags) 


(Stars) 


Source 


4118.593 


0.48 


1.873 ±0.018 


2.628 ± 0.036 


-0.055 ± 0.004 


SMARTS 


4121.529 


0.57 


1.888 ± 0.020 


2.571 ± 0.038 


-0.049 ± 0.004 


SMARTS 


4125.526 


0.66 


1.863 ±0.017 


2.602 ± 0.033 


-0.047 ± 0.004 


SMARTS 


4128.574 


0.51 


1.900 ± 0.024 


2.696 ± 0.052 


-0.047 ± 0.004 


SMARTS 


4133.581 


0.45 


1.858 ±0.017 


2.557 ±0.033 


-0.051 ± 0.004 


SMARTS 


4136.572 


0.49 


1.873 ±0.015 


2.578 ± 0.030 


-0.068 ± 0.004 


SMARTS 


4301.869 


1.05 


1.815 ±0.008 


2.565 ±0.012 


0.006 ± 0.003 


SMARTS 


4307.824 


0.70 


1.794 ± 0.009 


2.571 ±0.016 


-0.041 ± 0.003 


SMARTS 



Note. — HJD is the Heliocentric Julian Day - 2450000 days. The goodness of fit 
of the image, X 2 /Ndof, is used to rescale the formal uncertainties when greater than 
unity. The QSO A&B columns give the magnitudes of the quasar images relative to the 
comparison stars. The (Stars) column gives the mean magnitude of the standard stars 
for that epoch relative to their mean for all epochs. A few points in the lightcurves (in 
parentheses) were not used in the analysis. 



